Non-perturbative analysis of the Gribov-Zwanziger action 
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Abstract. In the non-perturbative regime the usual gauge fixing is not sufficient due to the Gribov problem. To deal with 
it one can restrict the integration in the path integral to the first Gribov region by using the Gribov-Zwanziger action. In 
its local form it features additional auxiliary fields which mix with the gluon at the two-point level. We present an explicit 
infrared analysis of this action. We show that from the two possible scaling solutions obtained previously only one remains: It 
coincides exactly with the results from the Faddeev-Popov action, i.e., the ghost propagator is infrared enhanced and the gluon 

propagator infrared suppressed and the corresponding power law behavior is described by only one parameter K = 0.5953 

This corroborates the argument by Zwanziger that for functional equations it suffices to take into account the appropriate 
boundary conditions and no explicit restriction in the path integral measure is required. 
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Understanding confinement, i.e., the absence of quarks 
and gluons from the physical spectrum, has been a long- 
standing challenge for the theory of strong interations. 
In the course of time several seemingly different con- 
finement scenarios have emerged. By now it has become 
clear that different explanations are not mutually exclu- 
sive but are only different manifestations. Especially in 
the last fifteen years functional methods have provided 
new insights into some of these scenarios, as the basic 
quantities in these approaches, the Green functions, give 
direct access to some of them, see, e.g., |1!|,|2|]. Because 
of its accessibility mostly the Landau gauge has been in- 
vestigated in this context, see, for example, 010,15100]. 

A particular gluon confinement picture is due to Gri- 
bov and Zwanziger [8, 9], where the gluon is removed 
from the physical spectrum due to a strong suppression 
of its propagator at low momenta. The argument is based 
on the insufficient gauge fixing implemented by the 
Faddeev-Popov procedure [8]. While this is not harm- 
ful for perturbation theory it is important for the non- 
perturbative regime. The so-called Gribov-Zwanziger ac- 
tion incorporates a refined restriction in field configura- 
tion space in an attempt to avoid this problem (fl |9[] : In- 
stead of integrating over all configurations fulfilling the 
Landau gauge condition dA = one takes into account 
only configurations from the first Gribov region £2 de- 
fined as {A\dA = 0, —dD > 0}. As a result already the 
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tree level gluon propagator is confined because it van- 
ishes at zero momentum: In this case positivity, a re- 
quirement for any physical particle, is violated maxi- 
mally [10]. The ghost propagator on the other hand is 
found to diverge more strongly than the tree-level part, 
namely like 1 / p 4 JH, [Till . This result is only obtained 
when enforcing the horizon condition, which arises as an 
additional condition in the gauge fixing procedure. Also 
in our approach it plays a pivotal role. 

Although the insufficiency of the standard gauge fix- 
ing has been long known, Dyson-Schwinger equations 
(DSEs) relied up to now on the Faddeev-Popov action. 
To understand why this is working we refer to a conjec- 
ture of Zwanziger Deriving DSEs with an inte- 
gration in field configuration space formally restricted to 
the Gribov region does not change their form. The rea- 
son is that cutting the integration at the Gribov horizon 
does not introduce any boundary terms as the Faddeev- 
Popov determinant, which appears in the integrand of the 
path integral, vanishes on the horizon. Hence one can use 
the Faddeev-Popov action but has to take into account a 
proper boundary condition. 

The original Gribov-Zwanziger action contains a non- 
local term 01, which can be localized with auxiliary 
fields 0}° h , 6)"b, and q>^ b 191]. The statistics obeyed 
by these fields is important: The first two are anti- 
commuting and the other two commuting. The additional 
adjoint color index accounts for the correct number of 
degrees of freedom. For our purposes the local action is 
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The anti-commuting fields r\" and rj" comprise the 
Faddeev-Popov ghosts, the real part of (p^ b and (pff 
and the fields (b^ b and d^f\ The index c runs from 

f (N 2 — 1) + 1. Treating all these fields with only one 
pair of fields is possible, because they all interact only 
via the Faddeev-Popov operator M = —dD. The massive 
parameter y is the Gribov parameter. Its value is fixed 
by the horizon condition. Although the part proportional 
to 7 breaks BRST invariance, see, e.g., ITllloi . this ac- 
tion is renormalizable. See 111411 and references therein 
for more details on the renormalization and the horizon 
condition. 

DSEs for this action can be derived with the Mathe- 
matica package DoDSE B15I1 which is of great help here, 
as due to the existence of mixed propagators even the 
two-point DSEs become quite complicated: The AA, VV, 
7777, VA and AV two-point function DSEs have 37, 5, 3, 
5 and 36 terms. The equations can be found in lfl6Tl . 

For the analysis of the equations we write the full two- 
point function := 8 2 T/8A8A as 



Yf v ac = 8" c p 2 ci(p 2 )P, v + 8 ac Ll{p 2 )p^p v . (3) 

For the two-point functions Fav '■= 8 2 T/8A8V and 
Tw '■= 8 2 T/8V8V we make ansatze which take only 
the tensors from the bare action into account: 
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The functions c;, i = {A,V,AV}, are non-perturbative 
dressing functions for the gluon, V field and mixed AV 
two-point functions. In principle more tensors are pos- 
sible. In order to obtain the corresponding propagators 
we have to invert the (2 x 2) two-point matrix T^^, 

ye{A,v}: 
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The propagator of the tj field is unaffected by the mixing: 



D 



r)r],ab ^p7]i),fl&^ — l %ab% 



cd 



= -8 ab 8 c 



■d 9 

p 



(9) 



For determining the behavior of Green functions at 
low momenta we study the DSEs in this limit. Then the 
integrand is dominated by small momenta and all dress- 
ing functions can be replaced by their infrared (IR) ex- 
pressions. For the dressing functions we make the an- 
satze Ci(p 2 ) — dj ■ (p 2 ) K ' and similarly for the vertices. 
The complete qualitative behavior is now described by 
the IR exponents jq. For the propagators we cannot say 
which part will be dominant in the IR, so we have to 
distinguish four scenarios. For every case one can then 
perform an analysis along the lines of ref. lfl7ll . where a 
general method for the calculation of possible IR scaling 
solutions is described. It was shown in ref. lfT6Tl that two 
scenarios are inconsistent, as the IR analyses yield con- 
flicting conditions for the IR exponents. Hence we con- 
cluded that only two possible solutions remain, which 
are, however, both in qualitative agreement with the re- 
sults from the Faddeev-Popov action and the IR leading 
diagrams are given by those which have a bare gluon- f]r\ 
or gluon-W vertex. Here we will present evidence that 
one of those solutions can also be discarded. 

This case was denoted III in IU6I1 and is given by 
Ka + Ky = 2k A v, i-e-, both parts of the determinant ap- 
pearing in the propagators, c^(p 2 )cv{p 2 ) +2Nc AV (p 2 ), 
scale equally. The scaling analysis resulted in K :— Ky — 
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1 , where 



K is a constant to be determined. For this one has to 
solve a system of equations consisting of the four two- 
point DSEs truncated to their IR leading diagrams. The 
unknowns are the combinations of the dressing func- 
tion constants I\ = dAd 2 ^, h = d^ydy, It, — dAd v and K. 
As the determinant remains unchanged here, the vari- 
ables appear in a more involved way than for the second 
case treated below, where the propagators have a simple 
dependence on the dressing functions. Consequently a 
complete analytic solution is not possible. However, one 
can solve for two variables analytically and try to find a 
solution numerically for the remaining equations. The re- 
quired analytic expressions can be obtained with the cur- 
rent version of DoDSE 01511 . h is then expressed analyti- 
cally from one of the equations as is I\ or 73. The results 



are plugged into the remaining two equations. Finally, 
they are plotted for the allowed values of the variables: 
K is known to be positive from the scaling analysis and 
it has to be smaller than 1 in order for the Fourier trans- 
formations of the propagators to exists. The values of I\ 
and I?, must be positive as the dressings dU and d-q are 
necessarily positive. It turns out that no solution exists. 
Hence this case is ruled out for the chosen ansatz. 

Therefore the only solution in our truncation is given 
by k := Kv = fCjj > 0, K A = -2k + j - 2 and Kav > 
—k+j—2 and the dominant part in the determinant is 
c a (p 2 ) c v {p 2 )- However, for the V propagator the second 
term is suppressed compared to the first one, so that its 
final form in the IR is 

n7v' abCd = 1 2-^S aC S bd g^ (10) 

As a consequence the dependence of the IR leading 
diagrams in the two-point DSEs of the V and A field 
on the mixed propagator vanishes. Even more, the equa- 
tions turn out to be completely the same in the IR as 
those obtained from the Faddeev-Popov action, because 
the contributions from the V field and the rj fields can- 
cel each other such that only contributions equivalent 
to those from the Faddeev-Popov ghosts remain. Hence 
also the result for K, which can be calculated analyt- 
ically, is the same: 0.5953... fl Etl- Furthermore, the 
IR exponents of all vertices are determined as K"2n,m = 
(n - m) K+ (1 - n) (f - 2) JH [H. With this informa- 
tion the IR exponent of the mixed two-point function can 
be calculated. It turns out that there are several possible 
solutions for K^v, the smallest one being 0.0668776. 

We want to stress that the solution obtained here solves 
the complete tower of DSEs in the IR and is the only 
possible solution. An important point is how the horizon 
condition enters. Without it we only could get the trivial 
result of canonical scaling of all Green function because 
of the bare two-point functions appearing on the right- 
hand sides of the two-point DSEs. However, solving the 
DSEs requires a renormalization condition which can be 
chosen as the value of the ghost dressing function at 
zero momentum. It was shown in ref. 10] that this value 
decides between the scaling solution (when it is chosen 
to be 0) and a family of decoupling solutions. Here the 
horizon condition explicitly forces us to use the former 
renormalization 0. [a [IH ■ This is how this additional 
condition enters and why we get the scaling solution 
and no decoupling type of solution like, for example, in 
H(l |2l|] . In principle the analysis also allows for other 
possibilities, e.g., for an IR divergent gluon propagator, 
but they are not realized. Discussions about the different 
solutions can be found, e.g., in iPyL I2HI22II23I1 . 

Our results are in qualitative agreement with those ob- 
tained in refs. l24l l25ll . where some propagators of the 



auxiliary fields are found to diverge like l/p 4 . This cor- 
responds to the value k = 1, which was already obtained 
by Gribov |Ht] in a perturbative analysis. The agree- 
ment we find in the IR with the scaling solution of the 
Faddeev-Popov action is another indication of the solid- 
ity of this confinement scenario, since it corroborates the 
conjecture of Zwanziger jj, ll2[] about cutting the inte- 
gration at the Gribov horizon. 
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